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Abstract 

Marcella has presented a straightforward technique employing the Dirac 
formalism to calculate single- and double-slit interference patterns. He claims 
that no reference is made to classical optics or scattering theory and that his 
method therefore provides a purely quantum mechanical description of these 
experiments. He also presents his calculation as if no approximations are em- 
ployed. We show that he implicitly makes the same approximations found 
in classical treatments of interference and that no new physics has been intro- 
duced. At the same time, some of the quantum mechanical arguments Marcella 
gives are, at best, misleading. 
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In a 2002 paper Marcella[l] pointed out that quantum mechanics textbooks typi- 
cally "finesse" discussion of the famous double-slit interference experiment by resort- 
ing to analogies from classical wave optics in order to obtain the standard results. He 
goes on to present a simple method of calculating a single- or double-slit interference 
pattern via the Dirac formalism and claims that this provides a quantum description 
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of the experiment, moreover one that should be accessible to upper- level undergrad- 
uates. It is true that the usual textbook discussions of the two-slit experiment are 
riddled with ambiguities and omissions [2]; and this is for the reason that a consistent, 
quantum-mechanical treatment of the two-slit experiment does not, to our knowledge, 
exist. Such a treatment would necessarily introduce an interaction between the inci- 
dent particle beam and slits that results in a probabilistic distribution of particles on 
the observation screen. The Hamiltonian governing this interaction would, in general, 
be quite complicated, and it is far from clear that even in principle it could resolve 
the "measurement problem," the abrupt change of a deterministic system to a prob- 
abilistic one at the instant of measurement. One might therefore question whether 
Marcella's approach, which includes no explicit interactions — in particular no solu- 
tion to a wave equation — can claim to give an adequate accounting of the two-slit 
experiment. Furthermore, at first glance he makes makes virtually no assumptions 
whatsoever and makes no reference to any approximations, leaving the impression 
that his result is exact. We now show that he has merely carried out the same clas- 
sical interference calculation that students encounter in first- or second-year physics, 
and with similar assumptions. 
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Figure 1: Marcella's setup for the one-and two-slit experiment. After passing through the slit(s), 
the particle is assumed to be scattered at an angle with a momentum p and a y-momentum of 
p y = ps'm8. Note that 6 is measured from the centerline. 
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We focus on the single-slit experiment; the double-slit goes through in exactly the 
same manner. The following is essentially Marcella's analysis, with his experimental 
arrangement shown in figure 1. Note that the angle 9 is drawn from the center of the 
single slit or the midpoint between the two slits. If the system, in this case a particle, 
is initially in a state then according to the Born postulate, the probability P that 
the particle is scattered with y-momentum p y = p sin 9 is given by 

P(Pv) = \<Pv\1>)\ 2 , (1) 

where the bracket represents the overlap integral between the initial and final states. 
In the position representation, the momentum eigenfunction for a free particle with 
momentum p y is 

<l>=(y\py) = ^e^. (2) 

The scattering amplitude (p y \ip) can be calculated by inserting a complete basis of 
position states between {p y \ and \ip) and integrating: 



oo 



l \VyW)=\ {Py\y)(y\i>)dy = — h- I e '*'ip(y)dy, (3) 



where in the Dirac formalism the position wave function ip(y) = {y\ip). To calculate 
P{p y ) from this expression we need to know ip(y) at the position of the slit(s). Taking 
the slit width to be a, Marcella makes the most natural assumptions: 

The first statement is merely that the intensity is assumed to be uniform over the 
slit, or that there is an equal likelihood of finding the particle anywhere within it. In 
other words, we take \ip(y)\ 2 dy = Cdy for constant C. Integrating over the slit gives 
the total probability of finding the particle P = Ca = I. Therefore C — 1/a and 
ip(y) = 1/y/a. 
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We now simply plug this expression for ip into Eq. (3) to calculate the scattering 
amplitude: 
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With the definition a = ap y /2h, Eq. (1) gives for the probability of detecting a 
particle at angle 9 

„. . a /sin a\ 2 
the standard result for single-slit diffraction. 

Upon being presented with this procedure, one's first reaction might well be, 
"Where's the physics?" The Dirac formalism is exactly that, a formalism, and does 
not in and of itself add any physics to the problem. We see that the above technique 
includes no wave equation to be satisfied and no interaction — in other words, no 
dynamics. Has a miracle occurred? 

In classical optics one generally treats the scalar diffraction problem by finding a 
solution to the Helmholtz equation: 

(V 2 + k 2 )ip = 0, (7) 

of which the time-independent Schrodinger equation is the special case k 2 = 2m(E — 
V)/h 2 , for constant E and V. The Helmholtz equation cannot be solved exactly 
except in nearly trivial circumstances, and so one must resort to numerous approxi- 
mations. The usual ones for scalar diffraction theory are: 

1) The "freshman physics approximation," or basic Huygens construction, which 
assumes that each point within the slit acts as a source of spherical waves ip(r) = 
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e tkr /r radiating only in the forward direction. Note that this wave function is a 
solution to the source-free Helmholtz equation; 



Figure 2: A source s emits a spherical wavefront whose unit normal e Q is at an angle 0' with 
respect to the x-axis. At the aperture the wavefront is diffracted to an observation point P. The 
normal to the aperture is n and the unit vector from P to a point on the aperture is e r . 

2) More sophisticated treatments include an "obliquity factor" in the Huygens 
construction; the solution is assumed to be of the form ip ~ f(9,9')e tkr /r, where 9 
and 9' are as shown in figure 2. To find the obliquity factor, one generally employs 
the Kirchoff approximation [3, 4, 5], which in the case of an aperture in a plane screen 
assumes that a) the wave function ip and dip/dn = Vip ■ n, for n the unit normal 
to the aperture in the direction of propagation, are zero everywhere except in the 
aperture; b) in the aperture ip and dip/dn take on the value of the incident wave. 
Any text on electromagnetism or optics, including those just cited, shows that with 
the Kirchoff approximation the value of ip at the observation point P is given by the 
integral 



where A is a constant and the integration is over the aperture. As one sees from 
figure 2, the Kirchoff integral can be written 







(9) 



where the obliquity factor is f(9,9') = (l/2)(cos# + cos#'). 
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It is well known that specifying both ip and dip/dn overdetermines the system and 
makes the Kirchoff approximation inconsistent; nevertheless it generally gives results 
that are in good agreement with experiment for small angles in the short wavelength 
limit. One can remove the inconsistencies by choosing one boundary condition or the 
other; that is, by specifying either ip within the aperture or dip/dn but not both. 
Thus we have the other two approximations: 

3) Dirichlet approximation, in which ip is taken to be the indicident wave within 
the aperture. For the case of a plane aperture, the obliquity factor becomes / = 
cos(0); 

4) Neumann approximation, in which dip/dn is specified within the aperture. Here 
for a plane screen the obliquity factor becomes / = cos(6>') 

None of these four methods lead to exact solutions of the Helmholtz equation and 
all differ in the limit of large scattering angles and long wavelengths. However, all 
four agree in the standard slit diffraction regime of normal incidence, small scattering 
angles and short wavelengths. In this limit, with the additional Fraunhofer approxi- 
mation that the distance to the screen is large compared to the aperture width, the 
Kirchoff integral for an infinite slit of width a in the y-direction becomes 

na/2 

^j(P) = C e- lksin9y dy, (10) 

J -a/2 

where now 9 is the fixed angle between P and the center of the slit. 

The question is, Has Marcella made implicit use of any of this formalism? At first 
glance it appears not; he does not employ the Schrodinger equation and makes no 
reference to any approximations. He is purely concerned with what happens to the 
particle at the slit and makes no statement about what takes place at the distant 
observation screen. This is similar to a "hidden-variable" point of view, that the 
particle is emitted from the screen with a definite momentum before a measurement 
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takes place. It is just this sort of interpretative statement that leads to analyses that 
violate Bell's inequalities and is certainly not the sort of concept that should be im- 
planted in the mind of a student of quantum mechanics. If, however, the discussion 
is to take place within the framework of standard quantum theory, we have the right 
to demand that the Schrodinger equation be satisfied, as well as quantum mechan- 
ics' standard postulates. Thus, Marcella must be implicitly solving the Helmholtz 
equation in some limit. 

We can uncover this limit by a close examination of Marcella's few assumptions. 
Note first that the assumption of equal intensity across the slit made in Eq. (4) is 
equivalent to the statement that the incoming wavefront is a plane wave. The normal 
vector n of a spherical wavefront impinging on a plane aperture would vary in angle 
with respect to the screen normal n across the aperture; hence the relationship of 
the area elements and intensity could not be constant on the screen. This does not 
necessarily imply normal incidence; however, since Marcella does not specify phases, 
it might as well be. That is, the intensity could either be \ip \ 2 from waves normal 
to the aperture, or |-0 cos(6' / ) | 2 = \ip \ 2 for non-normal wave function ip. But since 9' 
must be constant, it can be absorbed into the definition of tp. Furthermore, Marcella 
has assumed that tp — outside the aperture, but has not said anything about dip /On. 
Thus he has implicitly assumed a Dirichlet boundary condition, and so the obliquity 
factor should be cos(0), which has been set to one, i.e., the paraxial approximation 
has been used. 

At the same time, we notice that in the integration of Eq. (5), he has held 
p y = psin(#) constant. He justifies this by saying that he is calculating the scattering 
amplitude of a particle being emitted at the slit with a fixed angle 9. As mentioned, 
he makes no statement about what takes place at the distant observation point. In 
standard quantum mechanics, however, one would assume that the wave function 
could be emitted at different angles from the top and bottom slit or integrate over 9. 
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That Marcella uses a single angle (figure 1) drawn from the center of the slit, or from 
the midpoint of the double-slit configuration is permissible only in the Fraunhofer 
limit. Thus, at best, he is performing the integral in Eq. (10) 

However, one can see that his assumptions that ip — 1/ y/a and use of the plane- 
wave eigenstates in Eq. (5) is not consistent. All that is being done in Eq.(5) is to take 
the Fourier transform of ip(y) = 1/y/a. If we transform from position to momentum 
space and back again, we should end up where we started: ip'(y) = ip{y) = l/y/a. 
Thus, 

i pk m 

il/ (y) = -= / <P(k)e tky dk, (11) 
V 2tt J-k m 

where we have let the limits of integration be finite and where 

0( fc )_L f' 2 e -^y' d y'. (12) 

V a J -a/2 



Hence, 
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If k m — > oo, then the last factor becomes 2n5(y — y') and the result i])'{y) = 1/a/o as 
required. However, in Marcella's problem p y = hk y = psin(#) and there is a maximum 
value of k. In that case 

= i r/ 2 M ^ ,0) 
T^Va J -a/2 [y - y ) 

which is nonanalytic and must be approximated for a given value of y. For k m (y — 
y') « 1 one gets Marcella's result only when, interestingly, A = 2a, but in that case 
the assumption that k m (y — y') « 1 cannot always be met. 

In fact Marcella's approximations are no better than the freshman physics ap- 
proach, which ignores all boundary conditions. Consider figure 3, in which we divide 
up a slit of width a into 2N Huygens sources, N above the central source and iV 
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Figure 3: Consider a slit of width a and divide it into 27V sources of Huygens wavelets, N above the 
central source and N below. For parallel rays the phase shift with respect to the central source will 
be ±fc(a/2) sin9(n/N) for the n-th source. The phase shift at the end of the slit is A = ±k(a/2) sin# 

below. Assume that the screen is at an infinite distance so that the rays are parallel. 
At the observation point P on the screen the wave function is merely the addition of 
all the Huygens wavelets. Thus 

N 

^(P)=MP) e i(wt+ H (15) 

n=-N 

where the n th phase shift with respect to the central source is 5 n = k(a/2) sm(0)n/N. 
Write 5 n = a n sin9 with a n = k(a/2)n/N. Notice that a n plays the role of k or 
momentum, while 6 is the position variable. The spacing between adjacent n's is 
Aq = {ka/2N) = constant. Multiplying and dividing Eq. (15) by Aa gives 

^(P) = t^ e ianSine Aa, (16) 

n=~N 

where we have absorbed the maximum value (say) of the rapidly varying function e lult 
into the definition of ip . Absorb the first Aa into the definition of ipo as we U- Then, 
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in the limit of large N, 
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which is exactly the standard answer. We thus see that the freshman physics con- 
struction of adding together Huygens wavelets is really a Fourier transform, which 
is exactly what Marcella has introduced by chopping his wave function at the edge 
of the slits. Although he is going from position to momentum space, by writing the 
result Eq. (6) in terms of 0, we have the same position-space result just obtained. 

In sum, Marcella does make the valid point that quantum interference should be 
treated as a quantum phenomenon and quantum texts ought not immediately redirect 
the discussion to classical wave optics. But a more reasonable way to do this would be 
to simply show that the Schrodinger equation reduces to the Helmholtz equation, thus 
reducing the problem to one of classical scalar scattering with its concomitant approx- 
imations. This would also provide the opportunity of discussing relevant boundary 
conditions and to point out the difficulty of specifying them precisely in both the 
quantum and electromagnetic cases. As it stands, while Marcella's procedure is use- 
ful in giving students practice with the Dirac formalism, it has introduced no quantum 
physics into the problem other than setting p = hk, and has implicitly made all the 
assumptions that show this is indeed a problem of classical optics. That his result 
is the same as the one obtained by the simplest Huygens construction is merely a 
reflection of the fact that he has implicitly made the lowest-order approximations, 
where all methods converge to the same result. 
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